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We compute hadron masses and the lowest moments of unpolarized and polarized nucleon structure functions 
down to pion masses of 300 MeV, in an effort to make unambiguous predictions at the physical light quark mass. 



1. INTRODUCTION 

Understanding the structure of hadrons in 
terms of quark and gluon constituents (partons), 
in particular how quarks and gluons provide the 
binding and spin of the nucleon, is one of the out- 
standing problems in particle physics. 

Moments of parton distribution functions, such 
as (x) and gA = Am - Ad, are benchmark cal- 
culations in lattice QCD. To compare the lattice 
results with experiment, one must extrapolate the 
data from the lowest calculated quark mass to the 
physical value. Up to now, most results are for 
quark masses of about the strange quark mass 
and larger. A naive, linear extrapolation of (x) 
in the quark mass overestimates the experimen- 
tal number by « 40% [0. While most results 
are for quenched QCD, where one might expect 
that (x) is larger than the experimental value, re- 
cent unquenched results H,||,|[ indicate that this 
problem remains. 

It has been argued that a linear extrapolation 
must fail because it omits non-analytic structure 
associated with chiral symmetry breaking J5| . In- 
deed, chiral perturbation theory suggests a 
large deviation of (x) and <?a from linearity as the 



pion (quark) mass tends to zero: 



{x)NS = {x)ns( 
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* Poster at Lattice 2002, Cambridge, USA. 



where A is a phenomenological parameter. Equa- 
tion (1) fits both the lattice data and the exper- 
imental value ||. However, A is not yet deter- 
mined by the lattice data. To constrain this pa- 
rameter, and to perform an accurate extrapola- 
tion based solely on lattice results, data at smaller 
quark masses are crucial. 

2. THE SIMULATION 

To narrow the gap between the 'chiral' regime 
and the lowest calculated quark mass, we have 
started simulations at quark masses correspond- 
ing to about twice the physical pion mass. The 
calculations arc done for Wilson fcrmions in 
the quenched approximation. Improved Wilson 
fcrmions are known to suffer from exceptional 
configurations, which would forbid such a calcu- 
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lation at currently accessible couplings. 

We work at = 6.0. Our present data sample 
consists of: 
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The 24 3 32 and 16 3 32 lattices are from our pre- 
vious runs |7]]. On all our lattices m,L > 4, so 
that finite volume effects may be expected to be 
small. We did not see any at K = 0.1563. We are 
currently increasing our statistics at the smaller 
quark masses. 

3. HADRON MASSES 

We first looked at the chiral behavior of pion, p 
and nucleoli masses. We did not find compelling 
evidence for non-analytic behavior in any of the 
three cases. In Fig. 1 we plot the nucleon mass 
as a function of the pion mass. Quenched chiral 
perturbation theory (qCPT) predicts j§| 

m N = m° N + C 1 / 2 rrii T + G\ml + C 3 / 2 m 3 (3) 

with 

C 1/2 = -^{D-3F) 2 S»-0.5. (4) 

Empirically one finds, to a very high precision, 

m N = Jm° N 2 + Ci 774, (5) 

which is an analytic function in the quark mass. 
In Fig. 1 we fit (3) and (5) to the data. Both 
fits are hardly distinguishable in the range of the 
data. Fit (3) gives C 1/2 = 0.7(1) and X 2 = 0.8. 
Note that the coefficient C\/ 2 comes out posi- 
tive instead of negative as it should, according 
to qCPT. There is some preference for the single 
parameter analytic fit (5), which gives \ 2 = 0.4. 
What comes as a surprise is that the resulting 
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Figure 1 . The nucleon mass against , together 
with the phenomenological fit (5) (solid curve) 
and the chiral fit (3) (dashed curve). We have 
taken the nucleon mass at the physical pion mass 
to set the scale, using (3). 

physical nucleon masses (denoted by open circles 
in Fig. 1) still differ by more than 10 %. We con- 
clude that pion masses of 300 MeV are not small 
enough to sufficiently constrain the fit function, 
so as to give unambiguous results. 

4. STRUCTURE FUNCTIONS 

Let us now turn to the moments of non-singlet, 
unpolarizcd and polarized structure functions. 
We consider {x)ns first. On the lattice one does 
not compute (x) ns directly but jl[] 

R = (x)Nsm N , (6) 

which can be interpreted as the fraction of the 
proton's mass carried by the u quark minus that 
carried by the d quark. Because chiral perturba- 
tion theory does not provide much guidance at 
present quark masses, we seek a phenomenologi- 
cal extrapolation. In Fig. 2 we plot R as a func- 
tion of the pion mass. The data points have been 
properly renormalized and converted to renor- 
malization group invariant (RGI) numbers follow- 
ing gj. We find that the data lie precisely on a 
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Figure 2. The lowest moment of the non-singlet, 
unpolarized structure function against m w , to- 
gether with a linear extrapolation to the physical 
pion mass. Also shown is the experimental result 
(-#■) taken from 0. 

straight line, all the way from « 1 GeV down 
to our smallest pion mass at = 300 MeV. The 
slope comes out to be rj 0.4, which is close to the 
value 1/3 expected in the heavy quark limit. A 
linear extrapolation to the physical pion mass, 
R = R° + -R1/2 m n , gives a value for (x)^ 1 that 
is much closer to the experimental number than 
previous results. A fit of (1) to {x)ns still does 
not constrain A. We find A = 360^ MeV. Our 
present data cannot distinguish between the fit 
function (1) and (R° + R 1 / 2 m v ) / 'y 'm° N 2 + C\m\. 

Let us now consider gA, the axial vector cou- 
pling of the nucleon. This quantity describes 
the fraction of the proton's spin carried by the 
u quark minus that carried by the d quark. In 
Fig. 3 we plot gA as a function of the pion mass. 
Again, the data lie on a straight line, and we do 
not see any sign of non-analytic behavior as sug- 
gested by (2). 

5. CONCLUSION 

It appears that at pion masses of 300 MeV we 
are not yet sensitive to the predictions of chiral 



Figure 3. The axial vector coupling of the nucleon 
against m ff , together with a linear extrapolation 
and the experimental number (-fc). 

perturbation theory. Whether the situation will 
improve on finer lattices has to be seen. 

The numerical calculations were performed on the 
APE100 at NIC (Zeuthen) and on the Cray T3E 
at NIC (Julich). 
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